Weak force detection with superposed coherent states. 
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We investigate the utility of non classical states of simple harmonic oscillators, particularly a 
superposition of coherent states, for sensitive force detection. We find that like squeezed states a 
superposition of coherent states allows displacement measurements at the Heisenberg limit. En- 
tangling many superpositions of coherent states offers a significant advantage over a single mode 
superposition states with the same mean photon number. 
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I. INTRODUCTION 

Non classical states of light have received consider- 
able attention in the field of quantum and atom optics. 
Many non-classical states of light have been experimen- 
tally produced and characterised. These states include 
photon number states, squeezed states and certain en- 
tangled states. There are a number of suggested, and 
actual, applications of these states in quantum informa- 
tion processing including; quantum cryptography jl], 0L 
quantum teleportation|g, ||, p [ [?] , Bj, dense coding||]] 
and quantum communication]]! 0[ |ll| , ]12fl to name but a 
few. They have also been proposed for high precision 
measurements such as improving the sensitivity of Ram- 
sey fringe interferometry|^| and the detection of weak 
tidal forces due to gravitational radiation. In this paper 
we consider how non classical states of simple harmonic 
oscillators may be used to improve the detection sensi- 
tivity of weak classical forces. 

When a classical force, F(t), acts for a fixed time on 
a simple harmonic oscillator, with resonance frequency u> 
and mass m, it displaces the complex amplitude of the 
oscillator in phase space with the amplitude and phase 
of the displacement determined by the time dependence 
of the force j|3|. In an interaction picture rotating at 
the oscillator frequency, the action of the force is simply 
represented by the unitary displacement operator 



D(a) — exp(aa^ — a* a) 



(1) 



where a, a) are the annihilation and creation operators 
for the single mode of the electromagnetic field satis- 
fying [a, a'] = 1, and a is a complex amplitude which 
determines the average field amplitude, (a) = a. For 
simplicity we will assume that the force displaces the os- 
cillator in a phase space direction that is orthogonal to 
the coherent amplitude of the initial state, which we take 
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to be real with no loss of generality. The displacement is 



thus in the momentum quadrature, Y 



-i(a — ad). To 



detect the force we would need to measure this quadra- 
ture. If the oscillator begins in a coherent state |ao), (cko 
is real) the displacement D(ie) causes the coherent state 
to evolve to e lea °\ao + ie). The signal is then measured 
to be S — (Y out ) = 2e, while the variance in the signal is 



given by V = (O - (Y m 
ratio is hence 



= 1. The signal to noise 



SNR = 



S 



= 2e 



(2) 



which must be greater than unity to be resolved (the 
measured signal must be greater than the uncertainty 
of this quadrature in a coherent state). Thus we find a 
standard quantum limit for the weak force detection as 



e SQL 
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(3) 



II. 



WEAK FORCE DETECTION WITH 
SQUEEZED STATES. 



It is well known |T^| that this limit may be overcome 
if the oscillator is first prepared in a squeezed state (a 
uniquely quantum mechanical state) for which the un- 
certainty in the momentum quadrature is reduced below 
the coherent state level. For the case of an appropriately 
squeezed vacuum state 



2n) 



(4) 



where the mean photon number is given by 

n = A 2 /(l-A 2 ) (5) 

and A = tanhr (with r being the squeezing parameter). A 
weak force causes a displacement D(ie/2) on the squeezed 
vacuum. In this case the signal and variance for the mea- 
sured momentum quadrature is given by[[l6| 

S = {Y ou t) = 2e (6) 
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and hence a signal to noise ratio of SNR 
minimum detectable force is given by[fl6| 



e > 



1 

2~e r 



2ee r The 



(8) 



which for large squeezing corresponds to t m in > l/Ay/E. 
We see that squeezing provides an increased sensitivity 
that scales as \j\ffi. 

Following early work by Bollinger et al. [[l7) , Huelga et 
al. have shown that quantum entangled states can be 
used to improve the sensitivity of frequency estimation 
using Ramsey fringe interferometry. Can entanglement 
be used to improve the sensitivity for force detection? To 
begin, let us consider an entangled state of two harmonic 
oscillators, the two mode squeezed state, 



|V) = Vl-A 2 ^A>,rc) 



(9) 



71=0 



where \n, n) — \n)x®\n)2- The entanglement in this state 
can be seen in a variety of ways. Most obviously it is an 
eigenstate of the number difference operator a\a\ — a\a2, 
between the two modes, and in the limit of large squeez- 
ing, A — > 1, a near eigenstate of phase sum[Q. Alterna- 
tively we can consider the correlations between quadra- 
ture phase operators. In the limit of large squeezing 
(A — * 1), the state approaches a simultaneous eigenstate 
of both X\ — X2 and Y\ +Y2, which is the kind of state 
considered by Einstein Podolsky and Rosen|l9|]. This 
kind of correlation has been exploited by Furasawa et 
a/. ]20|] to realise an experimental teleportation protocol. 
With two oscillators, we need to specify how the weak 
force acts. We will specify that the force acts indepen- 
dently on each oscillator. To detect the force, consider a 
measurement of the joint physical quantity described by 
the operator Y\ + Y2. It is then straightforward to show 
that the signal and variance of the measured result, after 
the displacement, are given by 



S 



(Y, + Fo 



4e 



V = ((y 1 +f 2 ) 2 )-(F+^) 2 = 2e 



(10) 

(11) 



which gives a a signal to noise ratio of SNR = 2^/2ee r . 
The minimum detectable force is then e > l/(2\/2e r ) 
which is a y2 improvement over the single mode squeezed 
state. For large squeezing the minimum detectable force 
can be expressed in terms of the total mean photon num- 
ber for both modes. In this limit e m i n ~ l/(4^/n tot ). 
This is the same scaling as we found for a single mode 
squeezed state. The apparent improvement due to entan- 
glement is simply a reflection of the fact that we have a 
two mode resource with double the mean photon number. 

For the two mode squeezed state with the measure- 
ment scheme chosen, there is simple way to understand 
this result. The entangled two mode squeezed state, (||), is 
easily disentangled by the application of a unitary opera- 
tor of the form U = exp(— iTr(a\a2 + a±a\) / 4) , which does 



not change the total energy. We will refer to this unitary 
transformation as the beam splitter transformation as in 
the case that the two oscillator modes correspond to op- 
tical fields modes, this transformation describes the scat- 
tering matrix of an optical beam splitter. The resulting 
state becomes a (disentangled) product state of two sin- 
gle mode squeezed states (as in Eq.(|J)). The weak force 
now acts to displace each of the single mode squeezed 
states, each of which may be used to achieve the squeezed 
state limit for displacement detection. As there are two 
realisations of the measurement scheme there will be an 
additional 1 / V2 improvement in sensitivity simply from 
classical statistics. It is thus inaccurate to attribute the 
improved force sensitivity of a two mode squeezed state 
to entanglement when Y\ + Y2 measurements are per- 
formed. In assessing the limits to force detection using 
entangled states of N harmonic oscillators we thus need 
to consider if any apparent improvement could have been 
achieved simply by using N copies of an appropriate non 
classical state of a single harmonic oscillator. 

Of course it may not always be so obvious to transform 
an entangled state to a product of non classical states. 
Consider an entangled state of the form 



I*) = ^Tc n \n,n) 



(12) 



71=0 



This state is correlated in number, but unlike the two 
mode squeezed state, it is not necessarily a near eigen- 
state of phase sum. If we consider a measurement of 
Y\ + Y2 as previously, the signal and variance after the 
displacement are 



S = 4e 

V = 2 (1 + {<Ja + tfb) 



(a t 6 t +ab)) 



(13) 
(14) 



which gives an improvement in the signal to noise ratio 
when (aJ a +tf b) < (a^b^ + ab). A state like this, with cor- 
related photon number, is the pair-coherent (or "circle") 
state given by|^l|, p2[ 



I circle) 



\cte~ 



\ae 



(15) 



where |. . .)„ and |. . .)& represent coherent states in the 
modes a and b. Af is a normalisation coefficient and a 
the amplitude of the coherent state. This state can be 
written in the form (TL2|) with 



1 



(16) 



Here Iq is a zeroth order modified Bcssel function. This 
state cannot be separated into product states via beam 
splitter transformations. It is easily shown that the min- 
imum detectable force occurs when 



1 



+ n - 



a 



(17) 
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with the mean photon number being given by n = 
all (2a) / Iq (2a). A small improvement is seen for all 
a, with the minimum occurring at a = 0.85 (e m in = 
0.221108). As a — > oo we have t m in ~ * 0.25. In this 
optimal region the mean photon number is small. The 
measurement of Yi + Yq. is not optimal however because 
it is not a near eigenstate. 

It is likely that one can achieve a significantly bet- 
ter sensitivity if one changes their measurement quantity 
from Y\ + Y2 to a quantity that is a near eigenstate of 
(|l2|). For these correlated photon number systems this 
could require a measurement of the photon number dif- 
ference of ([l2]) which with current technology is quite 
unpractical. 



III. WEAK FORCE DETECTION WITH CAT 
STATES. 

Let us now turn our attention to a less straightfor- 
ward example. In the previous example two entangled 
harmonic modes, the two mode squeezed state, gave an 
improvement in the signal to noise ratio (compared to 
a single mode) of With an entangled state com- 

prised of more modes, an even better improvement may 
be achievable. However there is no simple way to gener- 
alise the two mode squeezed state to give an entangled 
state of many modes. We now consider another class of 
non classical states, based on a coherent superposition of 
coherent states, which can be entangled over N modes. 

Consider N harmonic oscillators prepared in the cat 
state 



IN 



where 



= JV + {\a,a, 



\a,a, ... ,a) 



n 



a)k 



-a)) (18) 



(19) 



is tensor product of coherent states and N is the normal- 
isation constant given by 



J\f - 



1 



v / 2 + 2e- 2Ar l Q l 2 



(20) 



We take a to be real for convenience. For a >> 1 this 
normalisation constant approaches 1/V2, and we hence- 
forward make this assumption. Parkins and Larsabal[p3| 
recently suggested how this highly entangled state might 
be formed in the context of cavity QED and quantised 
motion of a trapped atom or ion. 

To begin our consideration of these states, let us con- 
sider the case of a single oscillator (N = 1), 



1 



^/2 



a)) 



(21) 



where the mean photon number is given by n = |a| 
When a weak classical force acts on the state in Eq.p 7 



it is displaced by 



1 / 

cos 9 



(e- ,7m W \a + /3)+ e lImi - a ^\ -a + (3)) 
-«|-a» 



e" 7 a) 



1 sxa.il — ) 



(22) 



where 9 = —Im(a(3*) and we have defined the even (|+)) 
and odd parity (|— )) eigenstates 



|±) = - 7 =(|a)±|-a)) 



(23) 



Our problem is thus reduced to finding the optimal read- 
out for the rotation parameter 9 for a two dimensional 
sub-manifold of parity eigenstates. The rotation is de- 
scribed by the unitary transformation 



(24) 



where a x = |+)(— | + |— )(+| is a Pauli matrix. 

The objective is now to find an optimal measurement 
scheme to estimate the rotation parameter, 9, and thus 
the force parameter, e. The maximum sensitivity will 
occur when 9 = —Im(a/3*) is maximised for a given 
displacement. Having chosen a real, 9 is maximised by 
choosing (3 purely imaginary. This corresponds to a dis- 
placement D(/3) entirely in the momentum quadrature. 
Setting (3 = ie, we have 9 = ea. The theory of opti- 
mal parameter estimation p4| indicates that the limit on 
the precision with which the rotation parameter can be 
determined is 



{59f > 



1 



Var(a x )i 



(25) 



where Var{a x )i n is the variance in the generator of the 
rotation in the input state |+), which is simply unity. 
Thus we find that uncertainty on the force parameter is 
bounded below by 6e > l/(2a). It thus follows that the 
minimum detectable force is e m in > l/(2a), which may 
be written in terms of the total mean excitation number 
of the input state as 



e > 



1 



2y/E 



(26) 



where the mean photon number n = \a\ 2 . This measure- 
ment is at the Heisenberg limit. Comparison with the 
result for the single mode squeezed state shows a similar 
dependence on the mean excitation number however the 
squeezed state sensitivity is better by a factor 1/2. 
We can now consider a two mode entangled cat state. 



|-0>i =AA(k 



-a)) 



(27) 



However this state is easily disentangled with the unitary 
transformation 



7f T 

U(tt/2) = exp[-i~(a\a,2 



(28) 
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(for a quantum optical realisation this is a 50:50 beam- 
splitter) to produce the separable state 



\4>) 



i = 7Vi7V 2 (|a)i + | - a)i) ® (|a) 2 + | - a) 2 )(29) 

As in the case for squeezed states, we only need consider 
the force detection sensitivity for the state of a single 
oscillator. The minimum detectable force is given by 



e > 



1 



2V2n 



(30) 



Here we see the v2 improvement from classical averaging. 
For the N mode state given by eqn ( |l8| ) a linear trans- 
formation also exists to transform the N mode entangled 
state to a product state of single mode cat states. In this 
case the minimum detectable force using TV modes, each 
prepared in cat state with amplitude a, is 



> 



1 



(31) 



As each mode has a mean photon number given by n — 
a 2 , the total mean photon number used in the entire 
experiment is h to t = Na 2 , the minimum detectable force 
can be written as t m in > 1/s/fhot- We see from here that 
there is no real advantage in using entangled states with 
the measurement protocol outlined, as the improvement 
is only the standard statistical improvement that one gets 
from multiple copies of a single mode cat state produced 
by disentangling the state. 



IV. ENTANGLED CAT STATES 

A question that be asked is whether both entangle- 
ment and collective measurements allow one to increase 
the sensitivity of this displacement measurement past the 
limits shown above? To address this question let us con- 
sider again the N mode entangled cat state 

IV') = ~^(\ a > a ' ■ • ■ ' a ) + 1 ~ a ' _a ' ■ ■ • ' ( 32 ) 

where the total photon number of the entire state is 
ntot — Na 2 . The weak force acts simultaneously on all 
modes of this N party entangled cat state. It causes a 
displacement D(ie) on each mode in ( |32| ) resulting in the 
state 



AN 8 



V2 



\a + le, a + le, . . . , a + le) 



(33) 



-iN9 



V2 



a + ie, —a + ie, . . , , —a + ie) 



where 9 = ea. The theory of optimal parameter es- 
timation indicates that the limit on the precision with 



which the rotation parameter is given by (25) but where 



a x = a xi ■ The uncertainty in this force parameter 

is hence bounded by 



1 

Na 



1 



VNnTo 



(34) 



and is at the Heisenberg Limit. We observe a critically 
important extra y/N improvement due to the entangled 
state and collective measurement (projective measure- 
ments onto \a,a, . . . , a) — \ — a, —a, . . . , —a)) which can 
be seen over N individual copies of the state \a) + 1 — a), 
or a single mode state \y/n to t) + \ — \Jn to t)- For a large 
and finite ntot it seems optimal from that one should cre- 
ate highly entangled cat state with as many modes as 
possible while maintaining a»l. 

In our consideration so far we have not considered the 
effects of loss or decoherence on these highly non-classical 
states. Whether we are considering highly entangled cat 
states or large amplitude single mode cat states these 
are all extremely sensitive to small amounts of loss and 
decoherence. Error correction and avoidance techniques 
can be employed to reduce these effects but are beyond 
the scope of this paper. 



V. GENERALISED CAT STATES. 

In the example just discussed, maximum sensitivity re- 
quired the classical force to displace the cat states in a 
direction orthogonal to the phase of the superposed co- 
herent amplitudes. In general there is no way to arrange 
this before hand, as the phase of the displacement de- 
pends on an unknown time dependence of the classical 
force. However a simple generalisation of the previous 
cat states can be used to relax this constraint. Note that 
the cat states are parity eigenstates and are thus the 
conditional states resulting from a measurement of the 
number operator modulo 2, = a) a mod 2, on an in- 
put state | a) with a real. We are thus led to consider the 
conditional states for measurements of fix — a) a mod K . 
Such states have previously been considered by Schneider 
et al.pH, Given a result v = 0, 1, . . . , K — 1 for such a 
measurement, the conditional (unnormalised) states are 



K-x 



\K,u) = ex P 



At=0 



K 



\ae 



(35) 



which are eigenstates of e l27Taia / K with eigenvalues 

e -i2nv/K 

The case of K = 4 has recently been considered by 
Zurekpr| in the context of decoherence and quantum 
chaos. Assume that the oscillator is initially prepared 
in the state 



14,0) 



ia) + 



(36) 



with a real. Under the action of a weak force charac- 
terised by a complex amplitude displacement /3, the out- 
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put state is 



e l6 \a) 



e*+\ia) 



-i<t>\ 



(37) 



where 9 — aIm(/3) and <f> — oTte(/3). The state now 
carries information on both the real and imaginary com- 
ponents of the displacement due to the force which may 
be extracted by measuring the projection operator onto 
the initial state. In the limit that K » \a\ 2 >> 1, the 
initial conditional state is simply the vacuum state and 
we recover the usual standard quantum limit for force 
detection by number measurement |l5| . 



VI. DISCUSSION AND CONCLUSION. 

We now compare our results to the study of Ramsey 
fringe interferometry introduced by Bollinger et al. |lj 
and discussed by Huelga et al. . In Ramsey fringe in- 
terferometry the objective is to detect the relative phase 
difference between two superposed states, {|0), |1)}. that 
form a basis for a two dimensional Hilbert space. These 
states could be the ground and excited states of an elec- 
tronic dipole transition. The problem reduces to a quan- 
tum parameter estimation problem. The unitary trans- 
formation which induces a relative phase in the specified 
basis is U{6) = exp[i9Z] where Z = |1)(1|-|0)(0|. We are 
free to choose the input state \ip)i and the measurement 
we make on the output state, which is described by an 
appropriate positive operator valued measure (POVM). 

The theory of quantum parameter estimation^] indi- 
cates in this case that we should choose the input state as 
= (|0) + \1))/V2 and the optimal measurement is a 
projective measurement in the basis |±) = |0) ± |1). The 
probability to obtain the result + is P(+\0) = cos 2 9. In 
N repetitions of the measurement the uncertainty in the 
inferred parameter is 



N 



(38) 



which achieves the lower bound for quantum phase pa- 
rameter estimation. Repeating the measurement N times 
is equivalent to a single product POVM on the initial 

product state n»=i ®(|0)s + | l^i) / V 7 ^- However it was 
first noted by Bollinger et al.JlTj that a more effective 
way to use the N level systems is to first prepare them 
in the maximally entangled state. 



|^) = -1(|0) 1 |0) 2 ...|0) 



N 



IDilD 



|l)iv) (39) 



and subjecting the entire state to the unitary transfor- 
mation U(9) = Yli=i ex p(~i9Zi), the uncertainty in the 
parameter estimation then achieves the Heisenberg lower 
bound of 



1 

N 



Briefly let us instead consider N/2 maximally entangled 
pairs. In this case we can combine ( fic| ) at N = 2 with 
the square root statistical benefit of N/2 repetitions This 
yields 



= 1/2 x y / 2/N = 1/V2N indicating that pair- 
wise entanglement yields only a margined benefit com- 
pared to full N-wise entanglement of for the phase esti- 
mation. 

We will now show that the entangled state in Eq.(|39|) 
is in fact a cat state for a collective operator algebra. 
The Hilbert space of N two level systems is the tensor 
product space of dimension 2 N . The entangled state in 
Eq.(|39"|) however resides in a lower dimensional subspace 
of permutation symmetric states [£7|. These states con- 
stitute an N + 1 dimensional irreducible representation 
of SU(2) with infinitesimal generators defined by 



1 N 1 N 1 N 



(41) 



(40) 



where Z % = \\)^\\ - |0)<<0|, X t = |1>;(0| + |0) 4 (1|, Z t = 
i|l)i(0| - i|0)i(l|. The Casimir invariant is J 2 = ±(J 2 + 

Jy + J 2 ) with eigenvalue ^r(^r + !)• The operator J z 
has eigenvalues m = —N/2, —N/2 + 1, . . . , N/2 which is 
one half the difference between the number of zeros and 
ones in an eigenstate. It is more convenient to use of 
the eigenstates, |to)jv/2> of these commuting operators as 
basis states in the permutation symmetric subspace. In 
this notation the entangled state defined in Eq. ( |39| ) may 
be written 

W = 4=(l - N/2)n/2 + W/2) N/2 ) (42) 

In this form we can regard the state as an SU(2) 'cat 
state' for N two-level atoms. Hence it is straightforward 
to see that a single 2 N level atom can achieve the same 
frequency sensitivity. Their equivalence can be also be 
understood by noting that the sensitivity of such fre- 
quency measurements is proportional to the energy differ- 
ence of the states involved. What entanglement allows is 
for one to create an effective state without the need of re- 
sorting to create a superposition between certain ground 
state and a highly excited one. 

A closer atomic analogy to a single mode cat state 
would be a cat state for a single N level electronic system. 
For example we could consider the unnormalised state 
defined on a hyperfine manifold with quantum number 
F, \F)p + I — F)p. Such states have been considered in 
reference |28j . A similar state could also be generated for 
the large mag netic molecular systems considered in refer- 
ences |2^, ||CJ . The key is in how the resources can be 
distributed and what type of measurements one is trying 
to achieve. If the single molecule can only be prepared 
with a certain N, then an advantage can be gained for 
frequency measurements by entangling the state of many 
single molecule systems pi], |3^|. However if we restrict 
the total system to having a fixed N and we have enough 
control so as to be able to prepare the system either 



6 



as a single large SU(2) molecular state or many entan- 
gled smaller molecular states, then the same sensitivity 
is achieved for the high precision frequency measurement 
(this was not the case for weak force measurements). 

To conclude, we have in this article shown how su- 
perpositions of coherent states can be used to achieve 
extremely sensitive force detection. For a single mode 
state | a) + \ — a) we have found that the minimum de- 
tectable displacement for weak force measurements scales 
inversely proportional to the square root of the mean pho- 
ton number of the superposition of coherent states. This 
is the same scaling obtained by a single mode squeezed 
state and achieves the Heisenberg limit for single mode 
displacement measurements. 

What is potentially more interesting is that if we take 
a number of individual copies of a single mode cat states 
then we still achieve the inverse square root scaling with 
total mean photon number (hence effectively allow one 
to increase the mean number of particles). If one starts 
with an N mode entangled cat state, then simple lin- 
ear transformation can be used to turn this state into N 
copies of a single mode cat state and hence achieve the 
tmin ~ 1/ y/ntot sensitivity. This however is not optimal 
as it does not achieve the Heisenberg limit for multi- 
ple modes. To achieve this limit for weak force detec- 
tion one must use both an entangled N mode cat state 



and a joint collective measurement (between the various 
modes). Entanglement is the critical resource to achieve 
the best sensitivity for a fixed n to t- On the other hand 
we have shown for frequency (phase) measurements that 
the sensitivity previously offered by entangling N two 
levels atoms can be achieved with a single 2 N level atom. 
The key is that the sensitivity for the frequency measure- 
ment is proportional to the energy difference of the states 
involved and both the entangled resource and the super- 
position resource have the same energy difference. En- 
tanglement allows one to create an effective state without 
the need of resorting to create a superposition between 
certain the ground state and a highly excited one. 
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